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Trails (bond-avoiding walks) provide an alternative lattice model of polymers to self-avoiding 
walks, and adding self-interaction at multiply visited sites gives a model of polymer collapse. Re- 
cently, a two-dimensional model (triangular lattice) where doubly and triply visited sites are given 
different weights was shown to display a rich phase diagram with first and second order collapse 
separated by a multi-critical point. A kinetic growth process of trails (KGT) was conjectured to map 
precisely to this multi-critical point. Two types of low temperature phases, globule and crystal-like, 
were encountered. 

Here, we investigate the collapse properties of a similar extended model of interacting lattice trails 
on the simple cubic lattice with separate weights for doubly and triply visited sites. Again we find 
first and second order collapse transitions dependent on the relative sizes of the doubly and triply 
visited energies. However we find no evidence of a low temperature crystal-like phase with only the 
globular phase in existence. 

Intriguingly, when the ratio of the energies is precisely that which separates the first order from 
the second-order regions anomalous finite-sized scaling appears. At the finite size location of the 
rounded transition clear evidence exists for a first order transition that persists in the thermodynamic 
limit. This location moves as the length increases, with its limit apparently at the point that maps 
to a KGT. However, if one fixes the temperature to sit at exactly this KGT point then only a 
critical point can be deduced from the data. The resolution of this apparent contradiction lies in 
the breaking of crossover scaling and the difference in the shift and transition width (crossover) 
exponents. 

PACS numbers: 05.50.+q, 05.70.fh, 61.41.+e 



I. INTRODUCTION 

The canonical lattice model of the configurations of 
a polymer in solution has been the self-avoiding walk 
(SAW) where configurations are those lattice paths that 
could be generated by a random walk on a lattice that is 
not allowed to visit the same lattice site more than once. 
Considered as a static equilibrium statistical mechanical 
ensemble, self-avoiding walks display an excluded volume 
effect: they are swollen in size if compared to the set of 
unrestricted random walks of the same length. A com- 
mon way to model intra-polymer interactions is to assign 
an energy to each non-consecutive pair of monomers lying 
on neighbouring lattice sites. This prescription defines 
the interacting self-avoiding walk (f SAW) model, which is 
the standard lattice model of polymer collapse using self- 
avoiding walks. The collapse transition in TSAW models, 
the so-called 0-point of polymers, is a second order phase 
transition that has been well studied. The standard the- 
ory [IH3] of the collapse transition is based on the n — s- 
limit of the magnetic tri-critical tfi 4 — <fi 6 0(n) field the- 
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ory and related Edwards model with two and three body 
forces [H [5] , which predicts an upper critical dimension 
of three with subtle scaling behaviour in that dimension. 

A physically equivalent way of obtaining the excluded 
volume in a random walk model is to prevent the walk 
from visiting the same bond more than once. This 
weaker restriction leads to another class of lattice paths 
called self-avoiding trails (SAT). The interacting version 
of trails, customarily obtained by giving a weight to mul- 
tiply occupied sites, also presents a collapse transition. 
The literature contains various definitions [BHH] of single 
energy models of interacting self-avoiding trails (IS AT) , 
all weighting multiply visited sites in different ways (they 
differ in how sites visited more than twice are assigned 
the energy of interaction). Regardless, theoretical pre- 
diction [9] and the evidence [TUhT^] suggests that the 
collapse transition of the ISAT model is in a different 
universality class to that of ISAW, although there is no 
clear understanding of why this is the case if true. 

On the other hand, a two-dimensional model (triangu- 
lar lattice) of an extended ISAT (elSAT) model, where 
doubly and triply visited sites are given different en- 
ergies, say — E2 and — e 3 respectively, was recently [Ti] 
shown to display a rich phase diagram with first and sec- 
ond order collapse separated by a multi-critical point. 
The occurrence of the type of transition depended on the 
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ratio 



of the energies given to multiply visited sites of different 
degree. 

In conjunction with this study a stochastic process, 
known as kinetic growth trails (KGT), was also consid- 
ered, where configurations of trails are produced by a 
growth process. The static configurations produced are 
trails, but the trails of a fixed length are not all produced 
with equal probability. The equivalent static model is an 
elSAT with a particular value of k = k G — 4.15 at one 
particular temperature T — Tq. It was demonstrated 
[H] that this value of k = k G separates models where 
the collapse transition is first order (k > k G ) from mod- 
els where the collapse transition is second order (k < k G ). 
It was shown that for k < k G the second order transition 
was most likely in a single universality class, which was 
the same as the one in the (two-dimensional) ISAW 9- 
point collapse transition. Moreover, the KGT process 
mapped exactly onto the transition point separating the 
two lines of first and second order transition. Defining 
the temperature of the transition in the elSAT model to 
be T c (k), it was deduced that 

T c (k G ) = T G . (2) 

Also, importantly it was found that while the second or- 
der transition encountered for small k had a low tem- 
perature phase that was globular, the low temperature 
phase at large k was crystal-like, in that the trail filled 
the lattice. 

In this paper we investigate a counterpart to the tri- 
angular lattice elSAT model on the simple cubic lattice. 
One reason for choosing the simple cubic lattice is that 
the coordination number is the same as on the triangular 
lattice. This means that if one considers a kinetic growth 
process KGT on the simple cubic lattice it will map to 
precisely that same location in the phase diagram of the 
three-dimensional model as the corresponding triangu- 
lar lattice model described above. Also the KGT on the 
simple cubic lattice has previously been studied [T2]: it 
demonstrates a scaling behaviour consistent with it map- 
ping to a critical point in the phase diagram of elSAT. 
One might then speculate that the phase diagram of the 
elSAT on the cubic lattice has the same structure as on 
the triangular lattice. However, we find there are impor- 
tant differences, including anomalous finite size scaling 
behaviour around the KGT point, and no evidence of 
any crystal-like low temperature phase. 

The paper is set out as follows. We begin by recall- 
ing the theoretical framework of the collapse transition. 
We then move onto defining the models including the ex- 
tended ISAT (elSAT) model on a general regular lattice, 
specifically on the triangular and cubic lattices, the vari- 
ous canonical ISAT models and the KGT process. Before 
describing our results for the cubic lattice we summarise 
the findings on the triangular lattice. 



II. COLLAPSE TRANSITION AND 
CROSSOVER SCALING 

There are two related ways of describing the collapse. 
One is by understanding the change in finite length n 
scaling of key quantities such as the radius of gyration, 
or alternatively end-to-end distance, and partition func- 
tion as the temperature is lowered past the transition 
temperature T c . The second is the associated singular 
thermodynamic behaviour in the thermodynamic limit 
(n — > oo) of the free energy per step and in the internal 
energy and/or specific heat at that same temperature. 

Let us first describe the finite size scaling change. As 
the temperature is varied there is a collapse transition at 
some T — T c . For high temperatures (T > T c ) the ex- 
cluded volume interaction is the dominant effect, and the 
behaviour is universally the same as the non-interacting 
SAW problem. The mean squared end-to-end distance 
R\ and partition function Z n are therefore expected to 
scale as 

R 2 n ~An 2 " with v > 1/2, (3) 

respectively with estimates of v and 7 to be 0.5874(2) [7 
and 1.156957(9) [15] in three dimensions. 

When fixed at the transition temperature T = T c for 
ISAW in three dimensions, tri-critical field theory and the 
Edwards model [H [5] predict similar scaling forms with 
v = 1/2 and 7 = 1, though with additive logarithmic 
corrections. 

At low temperatures (T > T c ) it is accepted that the 
partition function is dominated by configurations that 
are internally dense. The partition function should then 
scale differently to that at high temperature, since a col- 
lapsed polymer should have a well defined surface (and 
associated surface energy). One expects [TB] in three di- 
mensions asymptotics of the form 

Rl~An^, (4) 

with fi s < 1. So the exponent v = 1/3 and the fractal 
dimension of the polymer becomes 3. 

This change in scaling behaviour is reflected in the 
thermodynamic limit. In the thermodynamic limit there 
is expected to be a singularity in the free energy, which 
can be seen in its second derivative (the specific heat). 
Denoting the (intensive) finite length specific heat per 
monomer by c„(T), the thermodynamic limit is given by 
the long length limit as 

C(T) = lim c„(T) . (5) 

n— y oo 

In general one expects that the singular part of the spe- 
cific heat behaves as 

C(T) ~ B\T C - T\~ a , (6) 
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where a < 1 for a second-order phase transition. The 
singular part of the thermodynamic limit internal energy 
behaves as 



U(T) ~ B\T C - T\ 



(7) 



if the transition is second-order, and there is a jump in 
the internal energy if the transition is first-order (an ef- 
fective value of a = 1). 

Moreover one expects crossover scaling forms [T7] to 
apply around this temperature, so that 

c n (T)^n a ^C([T-T c }n*), (8) 

with < <p < 1 if the transition is second-order and 



c n {T)~nC([T-T c ]n) 



(9) 



if the transition is first-order (that is, <p is effectively 1). 
Assuming standard crossover theory [IB] it was deduced 
in |17] that the exponents a and <f> are related via 

2-a=i (10) 

<P 

Regardless of whether the full crossover theory holds 
one can usually define an exponent ip from the shift of 
the transition temperature at finite length T c „ measured 
by, say, finding the position of a peak in the specific heat 
as follows 



Z n (u) 2 , . ..,Wq) 



and the probability distribution by 



p E ((t>n;U2, ...,w g ) = 



(11) 

Another exponent can be related to the width AT n of 
the transition. We shall define this exponent <J), via 

AT„ ~ fin"*. (12) 

We use the same symbol <p because in a crossover theory 
the exponent will be the same (f> as described in (J8|. For 
a first order transition one expects that both (p and ip are 
1. For a second order transition in the standard scaling 
picture one expects the two exponents to be equal, ip — <p, 
although we note a breaking down of the standard scaling 
has been already observed in [19) in high dimensions. 

For the ISAW model the tri-critical field theory expects 
the width (crossover) exponent <f> = 1/2 and shift expo- 
nent ip = 1/2 with logarithmic corrections present in the 
scaling forms because the system is at the upper criti- 
cal dimension. The value of the specific heat exponent 
a = is consistent with the scaling relation (10), as a ^ 
logarithmically divergent specific heat [H[S] is predicted, {Q) n = ~y~ < 5( < ^ > ™) 



III. MODELS 

A. The general extended ISAT model (elSAT) 

On a lattice with coordination number greater than 3 1 
the self-interaction in the trail model can be implemented 
in different ways, depending on how the weight associated 
with contact site depends on how many times the site 
has been visited. The canonical ISAT model, which has 
been defined differently by different authors, fixes the en- 
ergy associated with sites visited more than twice based 
upon the energy of doubly visited sites. The elSAT is a 
generalised model that allows for different (independent) 
energies to be associated with multiply visited sites of 
different multiplicities. 

Consider a regular lattice of coordination number 2q 
(5 £ N,g > 2) and the configurations <p n € fi ra of trails 
of length n (bonds) starting from a fixed origin. Let — ee 
be the energy associated with lattice sites that have been 
visited £ times by the trail. Now let mg, I — 1, ... q be 
the number of lattice sites that have been visited I times 
by the trail. Note that one always has J^e ^ m t = 71 + 
l 2 . Hence, to each of these contact sites is associated an 
Boltzmann weight u>g = e^ et with u)\ = 1, where (3 = 1/T 
is the inverse temperature in suitable units of the inverse 
Boltzmann constant. The partition function is then given 
by: 



m 2 (<j>n) 



,m q {<j> n ) 



m 2 {4>n) , m g (<t>n) 



(15) 



(16) 



Z n ( UJ 2, ■ ■ ■ 7 Uq) 

We define a reduced finite-size free energy per step as 



K n = - logZ n , 



(17) 



related to the usual free energy per step as —j3f n = &n- 
The average of any quantity Q over the ensemble of 
allowed paths <p n 6 Sl n of length n is given generically by 



• Ul 



m q (c/> n ) 



(18) 



</>„esi„ 



c n (T c )^C(lnn) 3 / 11 . (13) 

The prediction [H [S] for the shift is 

T^-T^Dn-WQnn)- 7 ' 11 . (14) 

On the other hand Shapir and Oono [9] have argued 
that the ISAT collapse transition should be tri-critical 
in nature, as self-avoiding walks. However they predict 
that ISAW and ISAT are in different universality classes. 
Importantly, while the upper critical dimension for ISAW 
is expected to be d u = 3, the Shapir-Oono field theory 
gives d u = 4 for ISAT. 



The thermodynamic limit in this problem is given by the 
limit n — > 00, so that the thermodynamic free energy per 
step foo is given by 



p/00 — K oo — lim k Ti 



(19) 



1 On a lattice with coordination number less than 4, a self-avoiding 
trail cannot visit the same site twice and therefore is just a self- 
avoiding walk 

2 We will not count the initial occupation of the origin as a visit. 
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This quantity determines the asymptotic behaviour of 
the partition function, i.e., Z n grows to leading order 
exponentially as /j, n with /i = cxp . 



B. Cubic and triangular lattice elSAT 

Both the triangular and simple cubic lattices have co- 
ordination number 6. Therefore only two weights, lv 2 and 
u>3, appear in the formulae (15 1 and (16). 



In order to explore the two-parameter space of the 
cISAT we define a one-parameter family of models with 
weights defined by: 



(tu 2 , o>3 ) = {u,uj k ) (fc-elSAT) 



(20) 



for any positive real value of the parameter k, that is, the 
energies obey 



e 3 = ke 2 



(21) 



in the fc-elSAT model. We set the energy £2 = 1 for 
convenience from now on. 

In this parametrisation we define a reduced internal 
energy per step and a reduced specific heat per step in 
the usual way via 



u n = dn n /d\oguj 



du n jd logo; 



(m 2 + km 3 ) r 



(22) 



((m 2 + km 3 ) 2 ) n - (m 2 + fcm 3 )^ 



(23) 



Let us define the collapse transition as occurring at 

lu = oj c {k) (24) 

so that we expect for fixed k that 

c n (uj)~n a *C((uj-uj c )n*), (25) 

if crossover scaling occurs. More generally the shift ex- 
ponent is defined by 



Dn 



(26) 



where oj Cj „ is the location of the peak of the specific heat, 
and the width exponent by 



Auj n ~ E n 



(27) 



where Aw„ is the width of the half-height of the peak of 
the specific heat. Both these exponents could depend on 
k. 



C. "Canonical" ISAT models 

The canonical model used by Doukas et al. P3] is one 
where every successive visit to a site adds an energy —e 



to the total for that site so that a g-times visited sites 
has energy — (q — l)e. Therefore the canonical model is 
defined by the weight parametrisation 



(cJ2,u> 3 ) = (^,<^ 2 ) 



(28) 



The canonical model corresponds to the case k — 2 in 
our family of interacting trails. 

A second "canonical" model used by Prellberg and 
Owczarek [5J [T2] has k = 3 so that 



(w 2 ,u> 3 ) = (w,oj 3 ) . 



(29) 



In fact it seems that Meirovitch et al. [6] used lo 3 = ui 2 , 
that is k = 1. Interestingly, all these models k — 1,2,3 
have been seen [Tl] to behave in the same way on the 
triangular lattice and, as we shall see, seem to behave 
in the same way on the simple cubic lattice (though the 
two- and three-dimensional models differ from each other 
in behaviour). 



D. 



Phase diagram of elSAT on the triangular 
lattice 



25/3 




FIG. 1: Phase diagram of elSAT on the triangular lattice. 
The dashed transition line for small u>2 (large k) separating 
the coil from the crystal-like phases represents a first-order 
transition and ends at the solid circle which is the kinetic 
growth point (KGT). This KGT point is conjectured to be 
multi-critical. The solid line separating the coil from the 
globule represents the #-like second order transition found for 
small k. The line separating the globule from the crystal-like 
phase is conjectured to be second order but not 6-like. 

The study by Doukas et al. [H] on the triangular lattice 
identified the kinetic growth point with a multi-critical 
collapse transition, being the meeting point of a swollen 
(coil), a collapsed and a crystal- like phase (see Fig.[T]). 

For small w 2 and 0,3 the trails present the usual swollen 
polymer phase where v = 3/4 (in two dimensions). For 
large enough w 2 , regardless of oj 3 , there is a collapse 
phase, as occurs in the ISAW model, and a transition 
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between the swollen and collapsed globule phases which 
seems to be 0-like, with <p — 3/7 as expected in two di- 
mensions. On the other hand, for large enough W3 the 
ensemble is dominated by crystal-like configurations that 
are space filling and internally contain only triply-visited 
sites. Between the swollen phase and the crystal-like 
phase the collapse transition is first-order. Separating 
this line of first-order transitions from the line of 0-like 
transitions is a multi-critical point. This point is pre- 
cisely the point (wj , w|) to which the kinetic growth pro- 
cess of trails maps. The crystal- like and globule phase 
is separated by a strong second order transition, much 
stronger than the 6*-collapse. 



IV. THE KINETIC GROWTH PROCESS (KGT) 

We now revisit the kinetic growth process [20] of trails 
(KGT). Consider a regular lattice of coordination num- 
ber 2q (q G N), and consider a stochastic process defined 
as follows: starting at the origin, a lattice path is built 
up step-by-step by choosing between available continuing 
steps from unoccupied lattice bonds with equal probabil- 
ity. This dynamic process produces lattice paths that are 
self-avoiding trails, moreover it is easy to show that, on a 
coordination 6 lattice, a trail <f> n of length n is generated 
with probability 



Pg( 



- - 

6 V 5 



n-l 



(30) 



This has to be compared with the probability distribution 
(16 1 of the equilibrium model with weights (5/3,25/3) 



PE 



5 25 
™ ; 3'T 



7(2. £2*1 
Zj n{ 3 , 3 ) 

from which we can deduce 



25 x 
3 



Pci^n) OC p E </>n 



5 25 
3' T 



(31) 
(32) 



Note that the normalisation is different since the sum 
over all walks of fixed length gives the probability of walks 
being still open in the case of the growth process, and 
unity in the case of the equilibrium model. 

We shall refer to the weight choice to which the KGT 
maps to as the Kinetic Growth point in the (0^2,^3) 
plane: 



(^2,^3) 



5 25 
3' T 



(KGT). 



(33) 



The Kinetic Growth point does not correspond to any 
point in any of the canonical parametrisations described 



above ( III C I but it belongs to our family of interacting 
trails ( 20 ) with k equal to 



log 25/3 
log 5/3 



4.15. 



(34) 



The KGT point has 



Tr 



1 



log 5/3 



1.957.. 



(35) 



In order to match the definition in ( 22 ) at the KGT point, 

u„ (04,0)3) = u* n (36) 
c n {u* 2 ,ul) = c* n , (37) 

we define analogues of the internal energy and specific 
heat by 



(m 2 + k G m 3 ) n 



(38) 



c , = ((m 2 + k G m 3 ) 2 ) n - (m 2 + k G m 3 ) n 
n 

Kinetic Growth trails have been subject of various 
studies both in low and high dimensionality and on a 
variety of lattices (see, for example, [TT1 IT21 [T4] ) . 

Simulations on the simple cubic lattice show a diver- 
gent specific heat with a logarithmic divergence different 
from the one predicted by the Edwards model. 

In particular it was seen that 



C*(logn) 



C 



with 



C = 1.0 ±0.5. 

whereas the Edwards model has £ = 3/11. 



(40) 
(41) 



V. RESULTS FOR THE CUBIC LATTICE EISAT 

A. The KGT point 

We begin our investigation by revisiting the KGT point 
on the simple cubic lattice directly by simulating the 
elSAT at that point. That is, we have simulated the 
KGT point of the elSAT model at 



(w 2 ,w 3 ) = (w^Wg) 



5 25 

3' y 



(42) 



We simulated 10 5 realisations of the kinetic growth 
process collecting 8.7 x 10 4 samples of length 2 20 . We 
then computed the energy u n and the specific heat c n . 
The results confirm what already was reported in |12) : 
u n behaves as 

u 

Un^Uao = , (43) 

y/n 

as indicated in Fig. [2j and the specific heat c„ diverges 
as a power of log n 



c (log n) 



(44) 



We calculated the local slopes of log c„ versus log log n, 
which are plotted versus 1 / log n in Fig. [3] We estimate 
C ~ 1.0 ± 0.5 as reported in |12) . 
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FIG. 2: Plot of the energy tt n of kinetic growth trails on the 
simple cubic lattice against 1/y/ri, along with a line of best 
fit. 
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10 5 iterations of generating about S n ~ 9 • 10 9 samples at 
length 10 3 . Following [21], we also measured the number 
of samples adjusted by the number of their independent 
growth steps, obtaining — 4-10 7 "effective" samples. 

FlatPERM outputs an estimate W n ^ m of the total 
weight of the walks of length n and fixed value of m. 
From the total weight one can access physical quanti- 
ties over a broad range of temperatures through a simple 
weighted average, e.g. 



E m ™ Wry 

Em Wn : m 



and 



J2 m W n ,r. 



(45) 

(46) 



The analysis of the scaling of the specific heat peak is 
done by calculating the location of the peak of the specific 
heat w Cj „ and thereby evaluating cf l eak = c n (w C! „). 



X X log c„ (KGW) 

— log cP eak (flatPERM) 

- - log cn{uj c ) (flatPERM) 
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FIG. 4: Comparison of the growth of the specific heat at 
the kinetic growth point (lower dashes and crosses) with the 
growth of the peak of the specific heat (upper dots) for elSAT 
with k = kc- One clearly sees different growth rates in the 
upper and lower curves, in particular the peak of the specific 
heat (upper dots) grows linearly on this logarithmic scale. 



FIG. 3: Estimating the exponent ( in Eqn. fl44| for KGT: Lo- 
cal slopes C, n of log c n versus log log n plotted against 1/ log n. 
Error bars are obtained by a simple re-binning procedure. 



The fc G -eISAT model 



We have first simulated the fc-cISAT model (20) with 
k = Icq (the correspondence with KGT occurs at w = 
5/3) using the flatPERM algorithm [2TII22] . We ran S = 



At the kinetic growth point, we confirm the slow 
growth of the specific heat seen in the previous section 
(c.f. Eqn. Q). The KGT simulations and flatPERM 
simulations evaluated at uj c coincide as shown in Fig. |4j 

Somewhat unexpectedly the height of the peak of the 
specific heat, c^ eafe , shows a different behaviour, diverg- 
ing as a power of n, as illustrated by the top curve in 
Fig. |1J In fact the divergence is linear, which indicates a 
first order transition. The nature of the finite-size transi- 
tion is therefore very different from the one inferred from 
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the kinetic growth process. 




0.25 



FIG. 5: Normalised energy density of states near the kinetic- 
growth point for elSAT with k = kc- The values of uj have 
been chosen such that both peaks have equal height. One 
clearly sees the build-up of a first-order transition with a well 
defined gap between the two peaks. This will result in a latent 
heat in the thermodynamic limit. 

We have therefore simulated this model at longer 
lengths (10 4 ) using a thermal implementation of flat- 
PERM at several fixed temperatures around the expected 
location of the transition. We have found that the en- 
ergy distribution shows indeed a double peak (see Fig. [5]) 
with the two peaks getting more and more definite as the 
length scale increases. 

We infer from this that a simulation at lj c does not 
"see" the second peak and hence only shows the observed 
much weaker divergence of the specific heat. 

This is consistent with a scenario in which the shift and 
width of the transition scale with different exponents as 
in Eqns. ((26]) and p7l 



Fig. [6] indicates that the peak position of the specific 
heat converges to the kinetic growth point uj c = 5/3 with 
an expected inverse square-root like scaling, while Fig. 
[3 shows that the width of the transition decreases much 
faster with an inverse linear scale. 
We thus conclude that 



Dn 



-1/2 



(47) 



where cj Cj „ is the location of the peak of the specific heat, 
and 



AuJr, 



En" 



(48) 



where Aw„ is the half-width (or, more precisely, width 
of the half-height) of the specific heat peak. 
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FIG. 6: The peak position of the specific heat w ,n versus 
1/y/ri for elSAT with k = kc- The location of the kinetic 
growth point at ui = 5/3 is marked with a large dot on the 
vertical axis. One can see that the peak position of the specific 
heat approaches this point in the thermodynamic limit. 
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FIG. 7: The half width the specific heat Auj n versus 1/n for 
elSAT with k = k G . 



Since the width decays much faster than the shift, one 
cannot see the first-order nature of the true thermody- 
namic transition in a kinetic growth simulation which is 
fixed at the transition point; it lies outside the crossover 
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region. This is graphically indicated in Fig. [8] 



Auj n 
< > 



FIG. 8: The observed scaling at k = kc is such that (u) c 



C. fc-eISAT with k> k G 

We next discuss the case of k > kc, where on the 
triangular lattice a first-order transition was found. 

We simulated k = 5, k = 7, and a model where only 
triply visited sites are weighted (this can be seen as the 
limit of letting k — > oo). We simulated trails with length 
up to 10 3 , collecting at that length S n ~ 7 ■ 10 8 , 10 9 
and 2 • 10 8 samples corresponding respectively to S^/^ ~ 
2 • 10 7 , 3 • 10 7 and 4 • 10 6 effective samples. 

The energy distribution displays a clear double peaked 
form (Fig. [9| which becomes more pronounced as k in- 
creases. Moreover, it sharpens as n increases, which is 
clear evidence of a first-order phase transition. We find 
that due to the initial build-up of the bi-modality the spe- 
cific heat actually seems to increase faster than linearly. 
The scaling of the shift and width of the transition here 
is no longer consistent with the scenario found at k = kc ■ 
We now find evidence for 1/n scaling of the shift of the 
transition, rather than the 1/y/ri scaling found at k = ka- 



li. fc-elSAT with k < k G 

Finally we focus on fc-elSAT for k < kg. We simulated 
the fc-elSAT model with k = 0, 1, 2, 3, collecting for each 
simulation order of S n ~ 10 8 samples (corresponding to 
S^ ff ~ 10 6 ) at the maximum length N = 10 3 . 

On the triangular lattice, for k < kc the collapse tran- 
sition is 0-like, and hence second order. A corresponding 
conjecture for the cubic lattice is the presence of a weak 
second-order transition with a logarithmically divergent 
specific heat c„ ~ (In np. 

The results for k = 1 and k = 2, which we report in 
Fig. |10[ are consistent with this prediction. However, we 
note that our estimated value £ = 2.25 ± 0.25 is outside 
the prediction £ = 3/11 for the 0-point. 

While models with k = and k = 3 show very strong 
scaling corrections which make the analysis inconclusive 
there, we expect the second-order phase transition sce- 
nario to extend to these values, and indeed to the whole 
range of values < k < kc- 




4 ■ 10" 3 




FIG. 9: Normalised energy density of states (k > fee) for 
k = 5 and 7, at the respective values of cj = 1.65 and 1.44, 
i.e. near the peak of the specific heat curves. 



E. Low temperatures 

Motivated by the results on the triangular lattice, we 
now investigate the possible presence of a crystal-like 
phase in three dimensions. 

In a crystal-like phase, in which the trail asymptoti- 
cally fills the lattice, the quantity 1 — 3it3, i.e. the pro- 
portion of steps that are not involved with triply- visited 
sites per unit length, should tend to zero as n — > oo. 
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Following the analysis in [2], we show in Fig. 11 the 
quantity 1 — 3ii3 against n~ 2 / 3 at two points in the pa- 
rameter region of the low-temperature phase. This is the 
expected order of finite-size correction due to the pres- 
ence of a surface in a compact low-temperature cluster. 
The parameters (a>2,W3) chosen are representative of re- 
gions in the (x>i , W3~plane for which we would expect to 
observe a significant difference. While there is a substan- 
tial difference in the asymptotic value of 1 — 3it3 in these 
regions, it is clear that 1 — 3«3 does not tend to zero in 
either. 

We have also looked for a signature of a transition be- 
tween low temperature phases and were unable to find 
any emerging transition, unlike for the triangular lattice 
(see Fig. 13 below). Of course we cannot exclude that 



1.50 1.55 1.60 1.65 1.70 1.75 1.80 1.85 1.90 1.95 

log log n 



such a transition may become apparent at longer trail 
lengths. 



VI. PHASE DIAGRAM 



FIG. 10: Specific heat peak scaling for k = 1 and k — 2: the 
specific heat diverges logarithmically with an exponent close 
to C = 2.25. 
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FIG. 11: 1 — 3u3, the proportion of steps not involved with 
triply- visited sites per unit length, plotted against n _2//3 , for 
{lo2,ujz) = (4,4) and (\/30,30), respectively. 



Based on this criterion, the investigation on the triangu- 
lar lattice identified two different regions of the collapsed 
phase: one in which 1 — 3tt3 tended to zero, and one in 
which it did not. The one in which it did was associated 
with the first-order transition for k > kc- 



To investigate further the full two-dimensional phase 
diagram, we ran two more sets of simulations. First we 
ran a two-parameter flatPERM simulation of the general 
cISAT model. A simulation of this type is limited by both 
time and memory requirement but we have been able to 
collect ~ 10 8 effective samples at length N = 512. 




FIG. 12: Density plot of the logarithm of the largest eigen- 
value X m ax of the matrix of second derivatives of the free 
energy with respect to u?2 and 0J3 at length n = 512. 

A density plot of the maximum fluctuations, calculated 
from the eigenvalues of the matrix of second derivatives 
of the free energy is shown in Fig. |12| Our inference for 



the finite size phase diagram is shown in Fig. 13 
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FIG. 13: Schematic of the observed phase diagram at length 
n — 1000. Dots indicate the location of peaks of fluctuations. 
Empty and filled dots indicate putative first and second order 
phase transitions, respectively, whereas small grey dots indi- 
cate crossover. The dot on the k = k — G-curve corresponds 
to the special transition described in Sec. |V A| 

Then we investigated the phase diagram along vertical 
and horizontal slices, that is to say with u>2 and W3 fixed, 
for n up to 1000. Observed phase transition points are 
included in Fig. [13] In the collapsed region we find evi- 
dence for a wide crossover region, but no evidence for an 
actual transition between two distinct phases. 

We have provided numerical evidence of that the two 
different types of transitions lead to a single collapsed 
globule-like phase at low temperatures. Our conjecture 
for the thermodynamic phase diagram is found in Fig. |14| 



VII. CONCLUSIONS 

We have investigated the collapse properties of an ex- 
tended family of interacting self-avoiding trails in three 
dimensions on the simple cubic lattice where doubly and 
triply visited sites are given weights oj 2 = e £2 / T and 

We have explored the general elSAT model by con- 
sidering a family of models satisfying W3 = ui\ with k 
positive real number. A kinetic growth process (KGT) 
of growing trails on the cubic lattice maps to one tem- 
perature of the kc — 4.15 equilibrium model. We find 
that the collapse is second-order if k < k G and first-order 
if k > kc- This resembles the triangular lattice finding 
(although the nature of the second order transition is dif- 
ferent). Interestingly, the low temperature phase for both 
k < kQ and k > kQ seems to be a disordered globular 
state. 



FIG. 14: Conjectured phase diagram for cubic lattice. Solid 
lines indicate phase transitions of the second order while the 
dashed indicates a first-order phase transition. The red dotted 
line only indicates the fc-elSAT model which passes through 
the KGT point. Empty and filled dots indicate an estimated 
location of the phase transitions coloured as in Fig. |13| 



Exactly at k — ka, the finite size scaling picture is 
particularly intriguing: the energy distribution displays 
a double peak form, indicating a first order type transi- 
tion but we observe different values for the shift exponent 
ip ~ 1/2 and the width (crossover) exponent ~ 1. The 
thermodynamic limit location of this first order transition 
when k = ka is the temperature Tq that maps to the ki- 
netic growth process. However, if one simulates directly 
at the point k — kc and T — Tq then the finite size scal- 
ing encountered is entirely second-order like and shows 
no sign of the first-order transition which dominates in 
the thermodynamic limit. This can be understood by 
appreciating that the finite size transition region shrinks 
quicker than its centre approaches the limiting tempera- 
ture. 

These results help to illuminate previous contradictory 
work for interacting trails on the diamond lattice [T2"ll23| . 
As the coordination number of the diamond lattice equals 
4, trails can only interact through doubly visited sites. 
The collapse point of interacting trails on the diamond 
lattice at u> = 3 was identified with the kinetic growth 
process. In it was shown that the scaling of the spe- 
cific heat at the kinetic growth points for the diamond 
and simple cubic lattices was indistinguishable. However, 
simulations of interacting trails on the diamond lattice 
showed the emergence of a first-order phase transition 
[23] ■ The scenario we describe here for KGT on the sim- 
ple cubic lattice clearly mimics these results. We are 
now able to understand the existence of both of these 
behaviours through the breaking of crossover scaling. 
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One last observation we can make is that interact- 
ing trails have been simulated in high dimensions El] 
and also demonstrate the breakdown of crossover scal- 
ing. The behaviour in high dimensions has been shown 
to be consistent with a self-consistent mean-field theory, 
which also displays bimodal energy distributions, though 
these do not lead to real first-order transitions in high 
dimensions. While that theory cannot not be applica- 
ble to k = kc in three dimensions (it predicts shift and 
width exponents both equal to 1/2), it would be inter- 
esting to formulate a mean-field theory of the transition 
that occurs for our elSAT model when k = feg. Conse- 
quently this may imply that the upper critical dimension 
for the k = kc cISAT models is less than three, and 
may in fact be two. Here we point out the numerical ob- 
servation of confluent logarithms in the two-dimensional 



kinetic growth trails [TTJ [T3] , lending further support to 
that assertion as logarithmic corrections typically appear 
at the upper critical dimension of a phase transition. 
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